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Abstract 

Simplified  models,  performance  correlations  and  rigorous  simulation  models  are  identified 
as  alternative  techniques  for  the  prediction  and  characterisation  of  the  performance  of  natural- 
circulation  solar  energy  water  heaters.  The  features  of  particular  examples  of  each  approach 
are  described.  In  their  different  appropriate  contexts,  each  technique  is  shown  to  provide  good 
agreement  with  measured  system  behaviour.  A  more  general  role  in  practical  system  design 
is  foreseen  for  rigorous  detailed  simulation.  Models  have  been  used  primarily  for  generic 
research  optimisations.  ©  2000  Elsevier  Science  Ltd.  All  rights  reserved. 


1.  Introduction 

The  heat  transfer  fluid  circulates  from  the  solar-energy  collector  to  the  hot-water 
store  in  a  thermosyphon  solar-energy  water  heater  by  buoyancy  forces.  Neither  a 
pump  nor  its  associated  control  sensors  and  actuators  are  required  to  circulate  the 
heat-transfer  fluid.  Such  solar  water  heaters  thus  have  both  lower  initial  and  running 
costs  compared  with  equivalent  forced-circulation  units.  Furthermore,  thermosyphon 
solar-energy  water  heaters  can  perform  as  effectively  as  forced-circulation  units  [4]. 

The  performance  of  natural-circulation  solar-energy  water  heaters  may  be  pre¬ 
dicted  analytically  by  three  broad  alternative  approaches: 
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Nomenclature 
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area  (m2) 

mean  specific  heat  capacity  (J  kg-1  K_1) 
diameter  of  the  collector  flow  channel  (m) 
collector’s  efficiency  factor 
solar  fraction 

total  daily  insolation  upon  collector  (Jm~2) 
height  (m) 

global  insolation  falling  upon  the  plane  of  collector  (Wm-2) 
length  of  flow  channel  (m) 
mass  of  water  (kg) 

mean  daily  thermo sy phonic  mass  flow  rate  (kg  s_1) 
number  of  times  the  water  contained  in  the  store  is  circulated 
through  the  collector  during  the  insolation  period 
integer 

mean  rate  of  delivery  of  useful  heat  from  the  collector  during  the 
insolation  period  (W) 
total  heat  delivery  (J) 

collector-plane’s  inclination  to  the  horizontal  (Degrees) 
temperature  (K) 
time  (s) 

steady-state  heat  loss  co-efficient  (Wm-2  K_1) 
dimensionless  groups  defined  by  Eqs.  (26),  (27)  and  (28), 
respectively 

absorptance  of  the  absorber  surface 
coefficient  defined  by  Eq.  (32) 

temperature  difference  between  the  fluid  at  the  inlet  and  outlet  of 

the  collector  (K) 

efficiency 

duration  of  insolation  period  (s) 
transmittance  of  collector  glazing 
T~Ta 


Subscripts 


a  of  the  ambient  environment 

c  for  the  collector 

d  refers  to  the  downcomer  pipe 

e  effective 

f  for  the  heat-transfer  fluid 

i  refers  to  collector  inlet 

L  refers  to  water  withdrawn  from  store 
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m  mean  value 

n  refers  to  number  of  passes 

p  for  the  collector  plate 

R  heat-removal  factor 

s  refers  to  the  system  as  a  whole 

T  for  the  hot-water  store 

u  refers  to  the  upriser  pipe 

z  for  the  volume  of  water  to  be  heated 

1  refers  to  the  fully-mixed  store,  immediately  prior  to  the  insolation 

period 

Superscripts 

'  denotes  absorber  plate  efficiency  factor 

denotes  average  value  of  the  variable 


•  simplified  models,  providing  average  or  cumulative  diurnal  predictions 

•  correlation  of  performance  characteristics  from  either  simulation  or  monitoring  of 
generic  systems,  and 

•  rigorous  dynamic  simulation  models  using  detailed  numerical  analysis,  providing 
hour-by-hour  or  more  frequent  performance  predictions 


2.  Simplified  models 

A  multiple-pass  solar  water  heater  is  defined  as  a  system  in  which  the  volume  of 
water  in  the  store  is  circulated  through  the  collector  many  times  each  day.  It  has 
been  observed  in  numerous  experimental  studies  that,  for  a  multiple  pass  thermosy¬ 
phon  solar-energy  water  heater,  the  difference  in  the  fluid  temperatures  at  the  inlet 
and  outlet  of  the  solar-energy  collector  is  constant  [5].  Whether  a  particular  system 
is  single  or  multiple  pass  depends  on  the  system’s  geometry  and  flow  resistance  as 
well  as  on  the  prevailing  insolation  level.  An  almost  constant  temperature  differential 
across  a  collector  (concomitant  with  multiple-pass  behaviour)  ensues,  provided  either 
a  critical  height  difference  between  the  collector  and  the  elevated  hot-water  store  for 
a  fixed  rate  of  heat  input,  or  a  critical  insolation  level  for  a  fixed  system  geometry 
has  been  exceeded  [6]. 

The  exhibition  of  an  almost  constant  temperature  differential  across  the  collector 
has  been  used  to  develop  a  particular  simplified  analysis  —  the  ‘A T  model’  —  for 
the  determination  of  daily  system  performance  [5].  The  advantage  of  this  particular 
approach  is  that  there  is  no  necessity  to  determine  the  frictional  flow  losses  within 
the  fluid  circuit  nor  to  perform  a  time-wise  analysis  of  the  temperature  variation 
around  the  thermosyphon  loop  to  ascertain  the  buoyancy  pressure  head. 
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To  develop  any  simplified  model,  many  assumptions  are  required.  For  the  ’A T 
model’,  these  are: 


1.  the  temperature  differential,  AT,  across  the  collector  remains  constant  throughout 
the  insolation  period.  The  duration  of  periods  at  the  start  and  end  of  each  insolation 
period  during  which  the  system  does  not  exhibit  this  behaviour  are  assumed  to 
be  much  shorter  than  the  insolation  period. 

2.  diurnal  mean  insolation  and  ambient  temperatures  can  be  employed  in  the  analysis 
without  introducing  significant  errors  in  the  calculated  values  of  the  net  daily 
heat  gain. 

3.  at  the  start  of  the  insolation  period,  the  hot- water  store  is  fully  mixed  at  a  tempera¬ 
ture  7V  During  the  period  of  insolation,  the  mean  store  temperature  rises  linearly 
with  time  to  a  value  Tz  at  sunset.  No  hot  water  is  withdrawn  from  the  system 
during  this  period.  The  latter  assumption  implies  that  the  system  is  intended  to 
produce  hot-water  for  withdrawal  either  that  evening  or  early  the  subsequent  day. 

4.  all  fluid  properties  are  independent  of  temperature. 


For  no  mixing  in  the  hot-water  store,  the  water  temperature  in  the  store  after  n  passes 
of  the  full  volume  of  the  store  through  the  collector  is: 

T„=Tt+nAT  (1) 

This  is  also  the  inlet  temperature,  Tf  i  for  the  (n+1)  pass.  If  the  total  number  of  passes 
is  N,  then  at  the  end  of  the  day 

T=T,+N  AT  (2) 


(T-ry 

2 


(3) 

(4) 


Alternatively,  when  the  store  is  fully  mixed,  the  rate  of  increase  of  the  store  tempera¬ 
ture  is: 


AT 

At 


thus 


m=Tx 


At  the  end  of  insolation  period, 
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The  inlet  temperature  at  the  collector  is  assumed  synonymous  with  the  tank  tempera¬ 
ture,  throughout  the  insolation  period.  The  diurnal  average  fluid  inlet  temperature  at 
the  collector  can  be  found  by  taking  the  diurnal  average  of  the  tank  temperature,  i.e., 


<J-TX) 

2 


(8) 

(9) 


Thus  the  mean  fluid  temperature  over  the  whole  day  at  the  inlet  to  the  collector  is 
given  by  Eq.  (4),  irrespective  of  the  degree  of  mixing  in  the  store.  The  useful  rate 
of  heat  provided  by  a  solar  energy  collector  can  be  described  by: 


Q = FrAc  [ha  -  U c(Tfj  -  Ta)] 


(10) 


where  FR,  the  heat  removal  factor  is  given  by 


mcf  I  i-UcAcF' 

^R=vrf(  1— exp 


AcUc 


me 


f 


(id 


The  rate  of  energy  delivery  by  the  collector  is: 


Q=rhcfAT 


(12) 


Substituting  from  Eqs.  (4),  (10)  and  (11)  into  Eq.  (12)  gives: 


r 


'X 


V 

m 


V. 


UCAT 


Iza-U 


Ti+E 

2 


(13) 


From  the  energy  balance  for  the  whole  system,  over  the  full  insolation  period  the 
average  collector  plate  temperature  is: 


Tp=Fr 


(14) 


for  which  it  is  assumed  that  the  heat  losses  from  the  up-riser  and  down-comer  are 
negligible,  thus: 
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T\\MzCf- 


U jA 


T  = 


2  / 


+mtSTCft 


UtAtQ 


TzUjAjO 


(15) 


m,Cf  -  2 


The  system’s  diurnal  heat-gain  efficiency  is  defined  as 
MZC/T—Tl) 


ns= 


IGA, 


(16) 


The  critical  minimum  vertical  height  of  the  down-comer  which  will  prevent  circu¬ 
lation  occurring  at  night  is  given  by  [5]. 


h,,= 


1 


-{(Td-rc)/csin  s+(Tz-Tu)hT 

U  u  1 d ) 

If,  for  the  particular  installation,  the  assumption  that,  at  night, 


rc=ra,  rd=ra+5  (in°C)  and  Tu= 


Tz+Ta 


(17) 


adequately  reflect  reality,  then  Eq.  (17)  becomes 


(T-T.J 

5  lc  sins+v  *  -hT 
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1 
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2 

(18) 


Using  initial  estimates  for  Tz  and  TP,  Uc  is  calculated.  Using  Eq.  (14)  a  new  value 
of  Tp  is  calculated,  this  process  is  repeated  until  convergence.  A  new  value  for  Tz 
is  found  from  Eq.  (15)  and,  using  the  convergence  value  of  Tp,  an  iterative  procedure 
is  used  to  obtain  a  solution  for  Tz.  If  mutually-incompatible  parameters  are  introduced 
into  the  procedure,  then  no  solution  will  be  found.  This  is  because  Eq.  (13)  cannot 
be  solved  as: 


r 


UcAT>Im-Uc< 


(Tx+Tz) 

2 


(19) 


When  compatible  values  for  the  overall  collector  heat  loss  co-efficient  are  chosen, 
the  condition  expressed  in  Eq.  (19)  will  not  arise;  for  example  if  the  value  of  Uc  is 
large  then  a  small  A T  would  be  expected. 

Theoretical  predictions  from  the  A T  model  have  been  compared  with  measured 
performance  data  from  two  identical  thermo sy phonic  solar-energy  water  heaters  [6]. 
Each  1.4  m2  collector  consisted  of  single-glazed  aluminium  flat  plate  absorbers, 
painted  matt  black:  their  fourteen  flow  channels  were  of  25  mm  internal  diameter. 
They  were  inclined  at  60°  to  the  horizontal  and  were  oriented  facing  south.  The 
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storage  tanks  were  each  1.6  m  tall  and  0.3  m  in  diameter,  giving  an  aspect  ratio  of 
5.3  and  capacity  of  0.1125  m3.  In  each  system  the  upriser  discharged  into  the  tank 
at  a  distance  of  30  mm  below  its  top,  and  the  downcomer  was  connected  at  the  level 
of  30  mm  above  the  tank  base.  The  data  acquired  included  the  time  histories  of 
temperatures  at  various  heights  in  the  collector  store  and  the  ambient  temperature 
using  copper/constantan  thermojunctions.  The  global  insolation  on  the  collector  sur¬ 
face  was  measured  with  a  Moll-Gorczynski  pyranometer,  mounted  in  the  same  plane 
as  the  collectors.  During  the  tests,  no  water  was  withdrawn  from  the  store.  Compari¬ 
sons  of  the  predicted  and  observed  values  for  the  mean  temperature  of  the  stored 
hot  water  are  shown  in  Fig.  1.  These  results,  obtained  for  a  wide  range  of  ambient 
and  initial  system  conditions,  show  good  agreement.  The  model  could  not  be  applied 
to  the  system  for  those  days  during  the  test  period  with  a  final  mean  store  temperature 
of  <38°C,  as  the  collector  temperature  differential  did  not  exhibit  a  constant  value 
throughout  the  majority  of  the  insolation  period. 


Experimental  data  with  bar 
indicating  possible  error  of  ±1°C 
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Fig.  1.  Comparison  of  the  experimentally-ob served  final  mean  temperature  of  the  stored  hot  water  with 
that  predicted  by  the  ‘AT”  simplified  model.  No  water  withdrawal  from  store. 
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3.  Correlation  of  performance  characteristics 

In  this  approach,  using  simplifying  assumptions,  an  analysis  of  the  thermal  behav¬ 
iour  is  undertaken  to  identify  dimensionless  groups.  The  developed  correlations 
between  these  characterising  groups  describe  long-term  system  performance.  As  the 
simplifying  assumptions  do  not  apply  in  all  cases,  the  functional  relationship  between 
the  dimensionless  correlating  parameters  obtained  analytically  is  not  that  observed 
in  reality.  The  actual  correlation  between  the  identified  grouped  parameters  is 
obtained  either  by  comparison  with  the  output  of  a  rigorous  numerical  simulation 
model  or  data  from  long-term  system  performance  monitoring. 

In  a  particular  approach  to  the  correlation  of  the  daily  performance  of  indirect 
thermosyphon  solar-energy  water  heaters  [2],  a  transient  heat  balance  was  undertaken 
on  a  generic  solar-energy  water  heater  (see  Fig.  2),  with  the  following  assumptions: 


-  Indicated  components 
of  energy  balance 


Fig.  2.  Holistic  heat  transfer  mechanisms  for  a  generic  thermosyphon  solar  energy  water  heater. 
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(i)  the  collector  and  storage  tank  both  have  the  same  mean  diurnal  temperatures, 

(ii)  insolation  and  ambient  temperature  are  assumed  to  remain  invariant  over  the 
day  at  their  respective  mean  measured  values, 

(iii)  the  total  draw-off  of  water  from  the  store  takes  place  as  a  single  event  at  the 
end  of  the  insolation  period, 

(iv)  water  is  drawn  off  at  the  mean  temperature  of  the  store  and  replenished  at 
the  mains  water  supply  temperature, 

(v)  tank  and  pipe  losses  are  negligible,  and 

(vi)  the  store  is  initially  at  the  mains  water  supply  temperature  at  the  onset  of  the 
insolation  period. 


An  instantaneous  heat  balance  on  the  system  gives,  with  assumption  (i) 


MrC^=^Ac[(m)ef  -  U(Tr-  Ta)] 


(20) 


Using  assumptions  (ii)  and  (v),  after  making  the  substitution  and  rearranging,  Eq. 
(20)  becomes, 


d0  ,  FRAcU<p  FRAc(za)cH 


MSCP6 


d  t  MsCp 

With  assumption  (vi)  as  a  boundary  condition,  the  solution  to  Eq.  (21)  is 

FRAc(za)eH-FRAcUe(Tm-Ta ) 


(21) 
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The  difference  in  the  mean  store  temperature  between  the  beginning  and  the  end  of 
the  insolation  period  is  therefore 


0-00  = 


F,W(Ta-)M-FKAM(Tm-Tj 
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From  assumptions  (iii)  and  (iv),  the  solar  fraction  can  be  expressed  as: 


/= 


MlCp(0-0,) 

Q 


(24) 


Substituting  Eq.  (23)  into  (24)  gives: 
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1 


(25) 


from  which  three  dimensionless  groups  can  be  defined  as: 
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For  the  particular  thermosyphon  solar-energy  water  heater  described,  inspection 
of  the  dimensionless  groups  X,  Y  and  Z  indicates  a  linear  relationship,  shown  graphi¬ 
cally  in  Fig.  3,  between  the  grouped  parameters  is  applicable  to  H(ra ),  FR ,  C/L,  Ac 
and  Mt  which  defined  the  system  configuration  and  thermal  characteristics.  The  per¬ 
formances  of  buoyancy-driven  systems  also  depend  on  other  factors  such  as  the 
position  of  the  store  relative  to  the  collector.  The  influence  of  these  other  factors 
may  be  encapsulated  in  the  relationship  between  the  store  and  collector  temperatures, 
that  is  the  validity,  or  otherwise  of  assumption  (i).  An  initial  investigation  indicated 
a  constant  temperature  difference  between  the  daily  mean  collector  (Tc)  and  store 
(TV)  temperature  over  the  insolation  period.  Fig.  4  shows  a  plot  of  the  daily  values 
for  Tc  against  Ts  for  the  month  of  August  indicating  a  linear  relationship.  The  best 
straight  line  fit  through  the  data  gave: 

T- 1.00  Tt=-  1.41  (31) 

with  a  correlation  coefficient  of  0.97.  Rearranging  Eq.  (31)  gives: 

7V-rc=1.41 
Thus,  more  generally 

Tt-T=P  (32) 

Introducing  into  the  analysis  such  a  constant  temperature  difference  /3,  between  the 
store  and  collector  results  in  the  dimensionless  group  Y  being  redefined  as, 


x 

(l-exp(-z)) 


Y 

and  —  correlated  over  a  period  of  250  days 


Fig.  3. 
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Fig.  4.  Correlation  of  the  mean  collector  and  hot  water  store  temperatures,  each  averaged  over  single 
periods  of  insolation. 
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Using  this  modified  value  of  Y  results  in  a  correlation  with  a  more  universal  appli¬ 
cation  potentially.  The  performance  of  any  system  can  then  be  determined  from 
several  daily  measurements  of  p.  Experimentally-determined  values  of  p  would 
enable  X ,  Z  and  a  modified  Y  parameter  to  be  used  to  predict  the  daily  performance 
of  the  system  under  the  applied  conditions  with  a  universal  correlation.  Alternatively, 
daily  values  of  p  could  be  generated  using  a  high-level  rigorous  simulation  model 
and  correlated  against  the  height  of  the  store  relative  to  that  of  the  collector.  The 
latter  approach  is  employed  by  [7]  to  develop  a  single  sizing-chart  nomogram. 


4.  Rigorous  dynamic  simulation  models 

In  a  rigorous  model,  the  coupled  energy  and  momentum  equations  which  describe 
the  transient  behaviour  of  a  natural-circulation  solar-energy  water  heater  are  solved 
numerically.  The  main  features  of  a  representative  example  of  a  dynamic  simulation 
model  [3]  are: 

(i)  a  finite  difference,  transient  heat  transfer  analysis  is  applied  to  the  liquid  circul- 
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ating  through  all  the  components  (i.e.,  collector,  upriser,  hot-water  store  and 
downcomer)  of  the  thermosyphon  loop, 

(ii)  the  thermal  capacitances  of  both  the  collector  plate  and  the  cover  are  included, 

(iii)  the  density,  specific  heat,  viscosity,  conductivity  and  Prandtl  number  of  the 
circulating  liquid  are  represented  as  second-order  polynomial  functions  of  tem¬ 
perature, 

(iv)  all  the  heat-transfer  coefficients  are  temperature  dependent  and  based  on  the 
ambient  and  mean  component  temperatures  and  are  updated  at  each  time  step  in 
the  numerical  calculation, 

(v)  the  store  model  includes  buoyancy-induced  mixing  between  stratified  layers 
which  occurs  when  warmer  fluid  is  introduced  below  a  cooler  layer, 

(vi)  a  transient  momentum  equation  is  employed, 

(vii)  friction  factors  are  calculated  using  correlations  appropriate  to  both  the  non- 
isothermal  thermally-destabilised  low  Reynolds  number  flow  and  isothermal 
developing  laminar  flow  in  straight  pipes, 

(viii)  empirically  determined  laminar  heat-loss  coefficients  for  the  pipe  bends 
are  included, 

(ix)  time  variations  of  insolation,  ambient  temperatures  and  hot-water  withdrawal 
encountered  by  an  installed  system  during  operation  [8],  were  used  as  inputs  to 
the  simulation, 

(x)  the  transmission  of  the  glass  collector  cover  is  determined  as  a  function  of 
the  sun  hour  angle,  and 

(xi)  a  two-dimensional  finite-difference  model  for  the  flat-plate  solar-energy  col¬ 
lector  was  employed.  The  main  advantages  of  this  over  previously  used  pro¬ 
cedures  are 

(a)  thermal  masses  and  temperature  profiles  for  each  of  the  fluid,  collector  plate 
and  glass  cover  can  be  described  individually;  and 

(b)  fewer  simplifying  assumptions  are  made  about  both  the  heat  transfer  processes 
within  the  collector  plate  and  between  the  collector  plate  and  the  fluid. 

The  equations  are  expressed  in  finite  difference  form  and  solved  simultaneously  using 
the  appropriate  boundary  conditions.  In  this  model  an  implicit  method  of  solution 
was  employed  in  which  N  simultaneous  equations  are  set  up  for  the  N  nodes,  and 
equations  were  solved  using  a  Gauss-Seidel  iterative  method.  For  this  method,  the 
solution  is  stable  unconditionally  and  the  size  of  the  time  step  is  only  limited  by  the 
required  accuracy  of  the  prediction  of  performance. 

Validation  was  undertaken  in  an  indoor  thermosyphonic  solar-energy  water  heater 
test  apparatus,  the  measured  global  insolation  was  simulated  experimentally  by  reg¬ 
ulating,  via  a  computer,  the  electrical  power  supplied  to  the  heater.  Similarly,  the 
thermal  load  on  the  system  was  reproduced  experimentally  by  the  controlled  opening 
and  closing  of  the  solenoid  valves  giving  a  draw-off  pattern  equivalent  to  that  pro¬ 
duced  by  the  occupants  of  the  house.  The  measured  insolation,  draw-off  pattern  and 
ambient  and  mains  temperatures  were  used  as  input  data  to  the  performance  simul¬ 
ation.  Measured  values  of  the  global  insolation  and  water  draw-off  patterns  were 
obtained  from  an  occupied  dwelling  in  the  south-east  of  England  to  which  a  thermo- 
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syphon  solar-energy  water-heater  had  been  installed  and  monitored  [7].  Fig.  5  shows 
the  variation  of  the  experimental  and  predicted  mean  store  temperatures  during  the 
day.  The  predicted  mean  store  temperatures  agree  to  within  0.7°C  of  the  correspond¬ 
ing  measured  values  over  the  period  of  time  from  00.00  hours  to  approximately  20. 
10  hours,  when  a  12  kg  water  draw-off  was  made.  After  this  draw-off,  the  predicted 
mean  store  temperature  fell  to  1.3°C  below  the  measured  value,  i.e.  the  predicted 
temperature  of  the  water  drawn  off  from  the  top  of  the  store  was  higher  than  the 
actual  temperature  of  the  draw-off  water. 

Good  agreement  is  shown,  as  illustrated  in  Fig.  6,  between  the  measured  and 
predicted  values  for  both  the  collector’s  inlet  and  outlet  temperatures  over  the  dur¬ 
ation  of  the  simulated  insolation  period.  Discrepancies  between  these  sets  of  values, 
particularly  at  the  end  of  the  insolation  period,  are  due  predominantly  to  the  uncer¬ 
tainty  in  predicting  the  frictional  flow  losses  at  the  very  low  Reynolds  number  flows 
encountered  during  reverse  circulation.  The  predicted  flow  rates  shown  in  Fig.  7 
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Fig.  5.  Diurnal  variations  of  experimental  and  predicted,  via  a  rigorous  simulation  model,  mean  store 
temperature  of  a  thermosyphon  solar  energy  water  heater. 
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Fig.  6.  Experimentally-observed  inlet  and  outlet  temperatures  of  a  flat-plate  collector  in  the  flow  circuit 
of  a  thermosyphon  solar-energy  water  heater  compared  with  those  predicted  from  a  rigorous  simulation 
model. 


correspond  to  within  7%  of  the  measured  values  over  the  insolation  period.  Flow- 
rates  were  not  measured  during  periods  of  reverse  circulation  as  no  signal  was  pro¬ 
vided  by  the  flowmeter  below  a  positive  flow  of  10~4  kg  s-1. 

The  large  fall-off  and  subsequent  recovery  in  the  measured  flow  rate  immediately 
after  water  had  been  drawn  off  from  the  store  was  predicted  correctly  via  the  math¬ 
ematical  model.  This  reduced  flow  rate  occurs  as  the  drawn-off  water  is  replenished 
by  cooler  water  (at  the  mains  supply  temperature)  introduced  at  the  base  of  the  store. 
This  cooler  fluid  enters  the  collector  and  the  resulting  increase  in  fluid  density  in 
the  risers  and  upriser  pipes  reduces  the  pressure  differential  between  the  fluid  in 
these  components  and  that  of  the  store.  A  subsequent  fall  in  the  flow  rate  ensues. 
The  reduced  flow  rate  caused  the  fluid  temperatures  in  the  collector  to  increase;  the 
fluid  density  falls  and  a  higher  rate  is  restored  rapidly.  The  total  thermal  output  due 
to  the  7  draw-offs  made  during  the  day  was  measured  as  4.27xl0~6  J.  The  corre- 
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Fig.  7.  Diurnal  variation  of  experimentally-ob served  thermosyphonic  mass  flow  rates  and  those  predicted 
from  a  rigorous  simulation  model. 


sponding  predicted  value  was  4.54xl0-6  J.  This  larger  heat  output  can  be  attributed 
directly  to  the  higher  degree  of  stratification  in  the  store  predicted  by  this  particular 
simulation  model.  Rigorous  dynamic  simulation  models  have  been  used  for  both 
direct  [2]  and  indirect  [9]  thermosyphon  solar  water  heaters,  the  latter  have  heat- 
exchangers  in  the  store  through  which  an  anti-freeze  solution  is  circulated  to  the 
collector.  Detailed  rigorous  simulation  models  have  been  employed  to  determine  the 
optimal  concentration  for  an  anti-freeze  solution  [8].  A  dynamic  simulation  model 
with  a  three-dimensional  representation  of  the  hot  water  store  has  been  employed  to 
analyse  the  effect  of  tank  geometry  on  thermal  stratification  [1]. 


5.  Conclusions 

The  three  generic  approaches  discussed  for  the  analysis  of  thermosyphon  solar 
energy  water  heaters  are  summarised  in  Table  1.  The  simplified  approach  described, 
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Table  1 

Salient  features  of  alternative  generic  approaches  to  the  analysis  of  thermosyphon  solar-energy  water  heat¬ 
ers 


Roles  Limitations 


Advantages 


Simplified 

models 


•  System  sizing  •  Only  apply  to  the  range  of  systems  and  •  Small  input  data 

operating  conditions  for  which  the  set 

simplifying  assumptions  are  valid. 

•  Often  require  experimentally-determined 
information  only  obtainable  after  the 
construction  of  the  unit. 


Correlation 

performance 

characteristics 


Long-term 

performance 

prediction 

Design  charts 

Economic 

optimisation 


•  Cannot  be  applied  reliably  to  those 
systems  for  whose  dimensions  and 
climatic  conditions,  a  correlation  has 
not  been  determined. 


•  Small  input  data 

set 

•  Wide-range  of 
systems  can  be 
considered 

•  Designer 
orientated 


Rigorous 

•  Engineering 

•  Require  increased  computing  power 

•  Accurate 

simulation 

optimisation 

•  Large  input  data  set 

transient 

models 

•  Establishing 

simulation 

long-term 

•  All  system 

performance 

configurations 

•  Determining 

can  be  included 

validity  of 

simplified 

models 

i.e.  the  A T  model,  is  very  quickly  executed.  Its  application  is  limited  by  the  need 
for  the  mean  diurnal  insolation  intensity  to  exceed  a,  usually  unknown,  threshold 
level  for  a  given  system  geometry.  The  principal  generic  shortcoming  of  such  simpli¬ 
fied  models  is  that  there  is  no  prior  way  of  determining  the  additional  performance 
parameter,  in  this  case  the  collector  temperature  differential  A T,  knowledge  of  which 
facilitates  the  simplification  of  the  analysis,  without  performing  either  an  experi¬ 
mental  measurement  or  a  rigorous  simulation.  However,  for  the  particular  case  of 
the  A T  model,  an  analysis  of  the  effects  of  changes  in  the  A T  value  employed 
revealed  that  though  the  fluid  mass  flow  rate  is  very  sensitive  to  this  factor,  the  final 
fluid  temperature  is  much  less  affected. 

The  performance  correlation  described  enables  the  long-term  monitored  perform¬ 
ance  of  an  indirectly  heated  thermosyphon  solar-energy  water-heater  to  be  related, 
on  a  daily  basis,  to  the  thermal  demands  made  on  the  system  and  the  prevalent 
meteorological  conditions.  The  correlation  exhibits  a  high  degree  of  linearity  over  a 
wide  range  of  climatic  conditions,  and  draw-off  patterns. 

Rigorous  dynamic  simulation  models  have  required  ‘expert’  users.  However  rapid 
and  continuing  advances  in  graphical  user  interfaces  and  pre-  and  post-processes 
together  with  the  ready  availability  of  computing  capacity  may  render  the  bespoke 
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design  of  large-scale  non-domestic  thermosyphon  solar  water  heaters  practicable 
using  such  models  as  tools. 
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